Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA.
نویسندگان
چکیده
In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA(TM) another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 A resolution crystal (structure the root-mean-square deviation of the 54 C(alpha) around the binding site and the inhibitor is 1.1 A). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.
منابع مشابه
Relative free energy of binding and binding mode calculations of HIV-1 RT inhibitors based on dock-MM-PB/GS.
Tetrahydroimidazo-[4,5,l-jk][1,4]-benzodiazepin-2-(1H)-one (TIBO) derivatives are important nonnucleoside human immunodeficiency virus-1 reverse transcriptase inhibitors (NNRTI). Several TIBO derivatives have shown high potency to inhibit reverse transcriptase (RT) and one (Tivirapine) has entered into clinical trials. The free energy of binding (FEB) is a numerical way to express the binding a...
متن کاملdMM-PBSA: A New HADDOCK Scoring Function for Protein-Peptide Docking
Molecular-docking programs coupled with suitable scoring functions are now established and very useful tools enabling computational chemists to rapidly screen large chemical databases and thereby to identify promising candidate compounds for further experimental processing. In a broader scenario, predicting binding affinity is one of the most critical and challenging components of computer-aide...
متن کاملIsothermal Titration Calorimetry and Molecular Dynamics Simulation Studies on the Binding of Indometacin with Human Serum Albumin
Human serum albumin (HSA) is the most abundant protein in the blood plasma. Drug binding to HSA is crucial to study the absorption, distribution, metabolism, efficiency and bioavailability of drug molecules. In this study, isothermal titration calorimetry and molecular dynamics simulation of HSA and its complex with indometacin (IM) were performed to investigate thermodynamics parameters and th...
متن کاملDesigning a new tetrapeptide to inhibit the BIR3 domain of the XIAP protein via molecular dynamics simulations
The XIAP protein is a member of apoptosis proteins family. The XIAP protein plays a central role in the inhibition of apoptosis and consists of three Baculoviral IAP Repeat domains. The BIR3 domain binds directly to the N-terminal of caspase-9 and therefore it inhibits apoptosis. N-terminal tetrapeptide region of SMAC protein can bind to BIR3, inhibit it and subsequently induce apoptosis. In th...
متن کاملg_mmpbsa - A GROMACS Tool for High-Throughput MM-PBSA Calculations
Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GR...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Journal of the American Chemical Society
دوره 123 22 شماره
صفحات -
تاریخ انتشار 2001